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Abstract. We show that nuclear motion of Rydberg atoms can be induced by 
resonant dipole-dipole interactions that trigger the energy transfer between two 
energetically close Rydberg states. How and if the atoms move depends on their 
initial arrangement as well as on the initial electronic excitation. Using a mixed 
quantum/classical propagation scheme we obtain the trajectories and kinetic energies 
of atoms, initially arranged in a regular chain and prepared in excitonic eigenstates. 
The influence of off-diagonal disorder on the motion of the atoms is examined and it 
is shown that irregularity in the arrangement of the atoms can lead to an acceleration 
of the nuclear dynamics. 
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1. Introduction 

Resonant energy transfer due to dipole-dipole interaction plays a crucial role in many 
areas of physics, chemistry, material sciences and biology Prominent examples are the 
light harvesting units of plants and bacteria pQ, aggregates consisting of organic dyes 
[21 [3] or molecular crystals Since the work of Franck and Teller [5] the theoretical 
framework to treat such processes is based on Frenkel's exciton concept [6]. However, 
internal and external degrees of freedom, e.g., internal vibrational modes of molecules 
and interactions with the surrounding solvent or proteins, have great influence on the 
properties of those systems and render a complete theoretical description challenging 
[7]. Hence, a prototypical system for a systematic study of the effect of dipole-dipole 
interaction should have (i) a simple structure, (ii) should not couple strongly to the 
environment and (ii) should be easy to manipulate, e.g. in terms of its interaction 
strength or geometrical arrangement. 

Rydberg atoms in ultracold gases are promising candidates to meet these 
requirements: First of all, their internal structure can often be reduced to a few relevant 
energy levels. Second, they posess long radiative lifetimes and their coupling to the 
environment via thermal collisions plays a minor role due to the low thermal velocities 
in ultracold gases. Finally, their mutual interactions can be controlled precisely, since 
the transition dipole moment of two adjacent Rydberg states increases quadratically 
with the principal quantum number. 

Consequently, clear evidences for resonant energy transfer due to dipole-dipole 
interaction have been observed in ultracold ("frozen" ) Rydberg gases [HI El [10, llj. There 
have been also several proposals how to study energy transfer in regular arrangements 
of Rydberg atoms p21 H3] • These geometries could in principle be produced by using 
optical lattices, or microlens arrays. Such one- dimensional chains have served as a 
fruitfull model to understand spectral and energy transfer properties of dipole-dipole 
interacting systems and to investigate the influence of disorder on their electronic 
transport properties. 

In studies on energy transfer the atoms/molecules are usually considered to remain 
at their initial positions. However, similar to the van der Waals interaction [14, 15J, 
the dipole-dipole interaction can lead to forces between Rydberg atoms, accelerating 
the particles. This effect has to be taken into account because it can lead to ionizing 
collisions and subsequent plasma formation, as demonstrated experimentally [THl CEO EE] ■ 
However, the induced forces might also be used to control the motion of the atoms, as 
recently suggested [T3] . 

In this work, we will investigate the motion induced by the resonant dipole-dipole 
interaction that triggers the energy transfer processes between atoms in energetically 
close Rydberg states, as observed in the studies on "frozen" Rydberg gases. To this 
end, we will assume that the Rydberg atoms are arranged in a linear array and focus on 
two distinct situations. First, we will consider a regular chain, in which all the Rydberg 
atoms are prepared such that they initially have equal nearest neighbor distance and 
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second, the situation when the Rydberg atom at one of the ends of the chain has a smaller 
distance to its neighbor than the atoms in the rest of the chain. The latter situation 
is a simple example, where the regularity of the mutual interactions in the system is 
broken. Such a system is often refered to as exhibiting "off-diagonal disorder". 

In both cases we will assume, that the system is initially prepared in one of 
the electronic eigenstates that form the exciton band. Experimentally, this could be 
realized by exciting atoms arranged in a microtrap array to a specific Rydberg state 
and subsequently promoting this system into an excitonic eigenstate of the chain by 
microwave excitation, similar to the procedure used in Refs. [T6l ITT] . 

A fully quantum mechanical treatment of this problem for long chains (i.e., many 
atoms) is computationally demanding. Therefore, we will use a mixed quantum/classical 
surface hopping method, in which the nuclear motion is treated classically while the 
electronic dynamics is treated quantum mechanically. 

The paper is organized as follows: In section [2] we give a brief description of the 
model system and the propagation scheme which we use. For the regular chain, we 
discuss analytical estimates for the initial forces as well as for the kinetic energy of the 
system and present numerical results in section [3] Having analyzed the basic properties 
of the atomic motion in the regular case, we proceed in section [4] with the case of broken 
regularity, by placing two atoms at one end of the chain closer together. Finally, section 
[5] presents our conclusions. 

Throughout this work atomic units will be used, unless stated otherwise. 



2. Rydberg atoms on a one dimensional chain 

2.1. Formulation of the model 

We consider a chain of N identical atoms with mass M and denote the position of the 
nth atom by R n . The distance R nm = \R m — R n \ between the atoms is assumed to be 
large enough that an overlap of their nuclear and electronic wavefunctions can safely 
be neglected. Since the focus of the present work is on resonant transitions between 
two energetically close Rydberg states, an essential state picture will be used, where for 
every atom only two electronic states are considered. We denote the energetically lower 
lying (Rydberg) state by \s) and the upper one by \p) and sometimes refer to the latter 
as the "excited state" or the "excitation". The energies of these states are e s and e p 
respectively, and u denotes the transition dipole between \s) and \p). By 

Kn) = \s---p---s) (1) 

we denote a state in which atom n is in state \p) and all others are in state \s). 
Without any interaction between the atoms the N states defined in ([I]) are degenerate 
eigenstates of the electronic system with energy e = e p + (N — l)e s which we choose to 
be the zero of energy and restrict our discussion to the space spanned by the states Q 
in the following. 
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The electronic energy transfer between two atoms n and m at distance R nm (Fig 
[TJa)) is mediated by the dipole-dipole interaction 

V nm {R nm ) , (2) 

where for the sake of simplicity we have ignored orientational effects. Due to this 
resonant dipole-dipole interaction the states given by Q are no eigenstates of the 
electronic system anymore. Moreover, as the coupling ^ depends explicitly on the 
positions of the atoms, their center of mass motion will also be affected by the dipole- 
dipole interaction, i.e., a force will be exerted on the particles that is proportional to 
the gradient of V nm with respect to the nuclear coordinates of the atoms. Consequently, 
both subsystems (electronic and center of mass) are coupled by ^ so that the motion 
of the particles will depend on the electronic state of the system and vice versa. 

The total Hamiltonian for a one dimensional chain of N Rydberg atoms is, thus, 
given by 

N n 2 

H(R U ...,R N ) = -52 —V\ n + H el (R 1 , ...,R N ), (3) 
where the electronic Hamiltonian, which depends on the nuclear coordinates, is 

N 

H C \R U ...,R N ) = - J2 V nm (R nm )\7r n )(n m \, (4) 

n,m=l 

and we set V nn = e = 0. 

2.2. Molecular dynamics with quantum transitions 

To study the motion of the Rydberg atoms induced by the dipole-dipole coupling we 
use a mixed quantum/classical method, where the nuclear motion is treated classically. 
To this end we denote the configuration of all atomic positions at a given time by 
p = {R\, . . . , Rn} and define adiabatic electronic states (p k (p) and energies U k (p) as 
the eigenstates and the corresponding eigenenergies of the electronic Hamiltoinan for a 
given p, 

H cl (p)Mp) = u k (p)Mp)- (5) 

Expanding the total wavefunction of the system in the adiabatic basis defined above, 

N 

*P(p,t) = J2 c k(t)Mp), (6) 

k=l 

and substituting ^ into the time-dependent Schrodinger equation of the full system we 
arrive at iV coupled evolution equations for the quantum amplitudes Ck, 

N 

ihc k = U k c k - ih^2 p ■ d kj Cj . (7) 
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Figure 1. (a) The p state is transferee! from atom to atom, due to the dipole- 
dipole interaction V nm . (b) Time evolution of the atomic distance in a dimer for 
different initial electronic configurations. The dotted lines are the trajectories, when 
the particles are propagated on one of the adiabatic potentials Uk with k = 1 and 
k = 2, respectively. The solid line is the trajectory, when the electronic excitation is 
initially completely localized on one of the atoms, the dashed curve, when 80 per cent 
of the electronic excitation is localized on one of the particles and both (cf. Eq{6]) 
have the same phase. Parameters are fj, = 1000 a.u., M — 11000 a.u. (i.e. Lithium). 



The terms that couple these equations are proportional to the velocities of the particles 
and to the so called "non-adiabatic coupling vector" dkj = { 4>k \ ^p(pj ) which, using the 
Hellmann-Feynman theorem, can be written as 

for j ^ k and dkk = 0. 

For a chain of length N = 2 ("dimer") all non-adiabatic coupling vectors vanish 
identically. In this case, the classical trajectories can be obtained by propagating the 
atoms on each of the two adiabatic potentials Uk according to the classical equations of 
motion 

MR n = F n = -V Rn ( <p k \H cl \ <Pk ) , (9) 

and averaging the results with respect to the initial probabilities \ck{t = 0)| 2 to be 
in one of the adiabatic states k. Figjl](b) shows the trajectories for a dimer, when the 
motion of the atoms starts on the adiabatic potential curves and two examples, in which 
the system is initially prepared in different superpositions of the electronic eiegnstates. 
From these simple examples it is already seen that the atomic motion induced by the 
dipole-dipole interaction depends crucially on the initial electronic state of the system. 

For chains with N > 3 the non-adiabatic couplings are generally nonzero. In 
this case we use a propagation scheme, in which an ensemble of trajectories is 
propagated, and each trajectory moves classically on a single adiabatic surface according 
to Eq. ([9]) except for the possibility of instantaneous switches among the adiabatic states 
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\4>k) due to the presence of the d k j- The algorithm used in this work to correctly 
apportion trajectories among the states according to their quantum probabilities with 
the minimum required number of quantum transitions is Tully's fewest switching 
algorithm [20]. For initial conditions that lead to fast collisions of atoms, the non- 
adiabatic couplings become significant, as they are proportitional to the velocities of 
the particles (c.f. Eq. [7]). This is manifest in large deviations of single trajectories from 
the ensemble average. Therefore, the mean trajectories for such initial conditions were 
determined up to the point t co n, at which the fastest collision event in a single realization 
occured. 

The numerical results that we will present in the following, have been obtained 
with the parameters M = 11000 a.u. (i.e., Lithium), /i = 1000 a.u. (corresponding 
to typical transition dipole moments between energetically close Rydberg states with 
principal quantum numbers n ~ 30, . . . ,40) and an initial nearest neighbor distance of 
x = 5 fim. 



3. Rydberg atoms arranged on a regular chain 

First, we will consider a chain consisting of N atoms, initially at rest with equal spacing 
R n ,n+i = x, and we will assume that the initial electronic state of the system is an 
eigenstate of H el . If one restricts the interaction of each particle to its nearest neighbors, 
it is possible to obtain a simple analytical estimate of the force ^ acting initially on a 
particular atom n, 

F n (t = 0) « 2 J^V(x) ( <p k |vr n > [( fc |vr n+1 ) - ( fe |7r n _i)] , (10) 

where we set ( 4> k \ttq) = ( 4> k \^n+i) = to account for the atoms at the ends of the chain. 
Thus, the initial force on particle n depends on the value of the electronic eigenstate 
(<f> k |7r n ±i) at both neighboring sites. 

For a regular chain with nearest neighbor interactions the eigenenergies U k and 
eigenstates | <p k ) are given by (see e.g. Ref. [TJ chapter 6) 

U k = -2V(x) cos ' ( U ) 

, , . / 2 / irkn \ , . 

\M=iN+i^ m [«rir'^ (12) 

The eigenenergies form a band, where for iV — > oo the ground state (k = 1) has energy 
—2V(x) and the highest eigenstate (k = N) has energy 2V(x). If the system is prepared 



in an electronic eigenstate \(f>k), the force on atom n can be obtained using (10) 



12 u 2 / nk \ ( 2nkn\ , , 

Hence, the force that initially acts on an atom decreases rapidly with increasing chain 



length N. This decrease of the magnitude is reasonable, since in the eigenstates (12) 



the excitation is delocalized more or less over the whole chain. Therefore, the excitation 
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is shared among all atoms in the chain, i.e., a particular atom carries only a fraction of 
approximately 1/N of the electronic excitation. 

The directions of the initial forces follow a simple scheme, that can be deduced using 
the symmetries of the system reflected in the nodal structure of the initial electronic 
wavefunction. Generally, the motion of the particles will be symmetric arround the 
center of mass coordinate of the chain (which is a constant of motion). In particular, 
for a chain with an odd number of atoms N the central atom (n = (N + l)/2) will 
always remain motionless, as the forces acting on it from both sides cancel each other at 
all times. Futhermore, there is one configuration, where all forces are directed towards 
the center of the chain and one, in which the accelerations are directed away from it. 
According to (13) this is the case, when the initial electronic wavefunction has no nodes 
at all (i.e., the state k — 1 at the lower edge of the band) and, when the sign of \4>k) 
changes from site to site (i.e. the state k = N with the highest energy). Generally, all 
forces reverse direction when changing the initial state from state k to (N+l) — k. A pair 
of atoms experiences an attractive force, if the wavefunction on the corresponding sites 
has equal sign and a repulsive force otherwise. The total acceleration of a particle due 
to its two neighbors is then determined by the relative magnitude of the wavefunction 
at neighboring sites (c.f. Eq. 10). 

For a chain with length N = 5, Fig. [2] shows the numerically calculated trajectories, 
for the system prepared in the five possible electronic eigenstates Uk, where now the 
interactions between all atoms have been taken into account. The graphs in the shaded 
areas show the corresponding eigenfunctions at the lattice sites ( (fik \n n ), where the site 
index n increases from bottom to top in each graph. 

The trajectories starting from the two energetically lowest eigenstates show 
relatively fast collisions eventually leading to ionization of the Rydberg atoms, whereas 
for the two highest eigenstates no collisions occur, although for k — 4 the trajectories 
of the three central atoms come close to each other. For k = 3 Eqjl3] predicts that no 
forces at all are present in the system, which seems to contradict the numerical result. 
However, the analytical estimate for the forces was derived under the assumption of 
nearest neighbor interactions only. In the general case of the "full" R~ 3 interaction, the 
excitonic band shows a marked asymmetry (Ui = —2.403 V(x) and Un = 1-803 V(x) for 
iV — > oo) and the phase of the eigenfunctions is shifted due to the additional interactions 
from more distant sites. Thus, the ( 0^ \iT n ) are nonzero at the sites next to the central 
atom, which in turn leads to a small but finite force between the particles. 

The dependence of the initial nuclear motion on the initial electronic state can also 
be seen in the time evolution of the kinetic energy of the system. For short times one 
can write E^ in (t) = (M/2) Z)^(t) with v n (t) « tF n /M. Inserting the initial forces (13) 
one finds 

. . 36 ( U 2 \ 2 9/ irk . t 2 , , s 

and E kin (t) ~ for k = (N + l)/2, if N is odd. We found good agreement of this 
formula with the numerical calculations shown in Fig. [3] for times t < 1 (is. The fast 
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Figure 2. Left column: Electronic eigenstates of a regular chain with N = 5 
atoms with increasing energies from top to bottom. Right column: Trajectories of 
the particles, when the system is initially prepared in the respective eigenstate on the 
left column. For details, see text. 



collsions in the trajectories for the two lowest initial electronic states lead to a rapid 
increase of the kinetic energy (Fig. km a)). This energy will eventually be available for 
ionizing the Rydberg atoms, when they approach each other. In contrast, E^ n increases 
much slower for the three eigenstates that do not show rapidly colliding trajectories 
(Fig. gb)). The slow "collision" that is seen in the trajectories for k = 4 manifests 
itself in a non-monotonic time evolution of Ey m with a minimum at the time, when the 
central atoms reach their minimal distance. 

4. Off-diagonal disorder - chain with "trap" at one end 

In this section we will study the effect, when the regularity of the chain is broken. Since, 
the regularity of the mutual atomic interactions V nm (for n ^ m) is destroyed, this is 



Motion of Rydberg atoms induced by resonant dipole-dipole interactions 



9 




t [|XS] t [JUS] 

Figure 3. Time evolution of the kinetic energy for a regular chain of N — 5 atoms, 
when the electronic state of the system is prepared in the five eigenstates. (a) Initial 
states k = 1 (solid) and k — 2 (dashed) that lead to fast collisions, (b) Initial states 
k = 3 (solid), k — 4 (dashed) and k — 5 (dash-dotted) that show a slow dynamics. 



an example of a system displaying off-diagonal disorder. We will focus on the specific 
case, where all atoms, but the last two have a spacing x and the distance between the 
last two particles (which we will call the "trap atoms") is y. 

In this case, there are generally no analytical solutions for the electronic 
eigenenergies and eigenstates, even if the interaction is restricted to nearest neighbor 
sites. However, if the eigenstates are determined numerically, Eq. 10 can still be used 
to estimate the forces acting initially on the particles. 

For a chain of iV = 5 atoms, Fig. [4] shows the trajectories together with the initial 
electronic eigenstates (<pk \^n), where we have initially set y/x = 3/5. In this case, the 
dipole-dipole interaction experienced by the trap atoms is approximately 5 times larger 
than by the particles in the rest of the chain. 

The eigenstate at the band edges {k = 1 and k = 5) are now strongly localized on 
the trap atoms and can approximately be written as ( | vr 4 ) ± | ir 5 ) )/ \/2, resembling an 
attractive and a repulsive dimer state, respectively. Consequently, the trajectories, when 
starting in these states, show, respectively, a strong attraction and a strong repulsion 
of the trap atoms. In the former case, the fast collsion is much better seen in the time 
evolution of the kinetic energy (Fig. [5^a)). 

Interestingly, if the system is prepared in the repulsive dimer state (k = 5), the 
"dimer character" is transferd from pair to pair. When two atoms come close to each 
other, they constitute a new repulsive dimer, such that the initial momentum of the first 
dimer is transfered through the chain. This effect can be enhanced, if the distance of 
the trap atoms is reduced further. Fig. [6] shows the time evolution of the kinetic energy 
for differnt ratios y/x ranging from 1 (i.e., completely regular chain) to 2/5, where the 
electronic eigenstate is partically completely localized on the trap. First, the kinetic 
energy of the system increases with decreasing distance of the trap atoms, as the initial 
electronic (potential) energy is dominated by the distance of the trap atoms. Second, 
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Figure 4. Same as Fig. [2] but for a chain with trap at the upper end. The ratio of 
the distance of the trap atoms to the nearest neighbor distance in the regular part of 
the chain is 3/5. 



.Ekin shows an increasingly non-monotonic behavior leading to pronounced oscillations 
for small y/x, where the minima correspond to times, when the distance of two atoms 
becomes minimal, i.e., a new dimer is formed. 

In the cases, when the system is prepared in on of the eigenstates corresponding 
to energies in the center of the exciton band, the excitation is delocalized over the 
regular part of the chain. The trajectories (Fig. [1]) and corresponding kinetic energies 
(Fig. [5]^a), (b)) show a similar behavior like for a completely regular chain, with the 
exception that in the slow collision process for k — 4 the Rydberg atoms approach each 
other so closely, that the assumption of non overlapping wavefunctions breaks down, so 
that this process has to be treated completely quantum mechanically. 
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Figure 6. Time evolution of the kinetic energy for a chain of N = 5 atoms with trap 
at one end starting from the k — 5 electronic eiegnstate for different ratios 7 = y/x of 
the distance of the trap atoms to the nearest neighbor distance in the regular part of 
the chain, (a) 7 = 1 (regular chain), (b) 7 = 4/5, (c) 7 = 3/5 and (d) 7 = 2/5. 
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5. Conclusions 

We have investigated the nuclear motion in a chain of Rydberg atoms, which is induced 
by resonant dipole-dipole interactions that trigger the energy transfer between two 
energetically close Rydberg states. 

The motion of the atoms depends crucially on the electronic state, in which 
the system is prepared. We have considered the particular case, when the system 
is prepaered in eigenstates of the exciton band. For a completely regular chain, we 
have given simple estimates for the initial forces and the time evolution of the kinetic 
energy for short times. In particular, we have shown that the magnitude of the forces 
acting initially is inversely proportional to the number of atoms, due to derealization 
of the electronic excitation over the whole chain. This derealization of energy is 
also responsible for phenomena like exchange narrowing [21] or superradiance found 
in molecular aggregates. 

Furthermore, we have studied the influence of so called "off-diagonal disorder" in 
a simple example, in which we have distrubed the regularity of the system by placing 
two atoms closer together than the particles in the rest of the chain. In this case, some 
electronic eigenfunctions are strongly localized on the perturbed sites and the degree of 
localization increases with increasing irregularity. In our example, where the irregularity 
leads to larger initial potential energy a much faster nuclear motion results in the system 
and a transfer of the initial momentum of the atoms forming dimer through the chain. 
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